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This study is targeted to the NLO corrections of multileg processes, very important for the LHC. 
Starting from the construction of Feynman diagrams, the analytical reduction of general one-loop 
integrals to scalar master ones, the calculation of color structures, manipulation of spinor lines 
and other amplitude constituents and finally phase space point selection are obtained by use of a 
program producing Fortran code for numerical calculation of one-loop corrections for processes 
like gg -> ttgg. 
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1. Introduction 

The era of LHC has finally arrived, after many years of waiting and anticipating, and all the 
particle physics community is waiting for the first results. There has been a lot of activity the last 
two years in the branch of NLO calculation. With almost all interesting cases of 2 — )• 3 processes 
having been calculated, the interest has moved towards 2^4. Processes Uke ;?/?—)■ tt + Z7Z> [|T], §, |3|] , 
V + 3 jet production [QISl, ^, ^, g^ — > bbbb ^ and most recently the case of pp ^ tt + 2 jets [^] 
have been studied in the last three years, providing the background for future discoveries. This 
study is also towards this direction. 

2. The method 



Our method is based on the Feynman diagrammatic approach. Which means that we calculate 
the total amplitude by taking into consideration the contributions of each of the one-loop Feynman 
diagrams that exist for the NLO calculation. The final program calculating the amplitude is in For- 
tran, although a variety of tools are used in combination (see Figure |2|): Diana [10| for constructing 
the diagrams. Form [|ri|] for manipulation of the Diana output (the steps of calculation in Form are 
colored red). Maple for optimization of the produced Fortran code (green color) and finally Fortran 
for the numerical evaluation (blue color). In the next session we will explain in more detail the 
complete calculation, which ends by providing the amplitude from an initial phase space point. 



2.1 The steps of the calculation 

As mentioned above, the first step is implemented by using Diana. We provide the program 
with the appropriate Feynman rules and ask for the construction of diagrams. For the case of 
gg — ;■ ttgg the total number of diagrams is 4510. Diana provides us with an output in Form which 
we further manipulate. The next step is to strip the one-loop amplitude from the color structures. 
In a few words using Form to manipulate the output of Diana we strip off and manipulate the 
color structures using color Algebra. We translate the Form result to Fortran code and end up 
with a Fortran array containing all the possible color information. For the case of gg — ;■ ttgg there 
are 50 color structures, which can be further divided in 4 main color groups. Representatives of 
each group are shown in Table [l[ All other color structures are produced from these ones through 
permutations of the external gluons. 



Color Structures 


Number of cases 


(r«ir«2r«3r«4)/,/2 


24 


j^|jaija2}(7«37«4)^^^^ 


12 


rr{r«ir«2r«3}(r«*)/,/2 


8 


rr{r«ir«2r«3r«^}5/,/2 


6 



Table 1: Color groups for the gg -^ ttgg case. The number of cases refers to all the possible permutations 
of the gluons 



The next step is the manipulation of the color stripped amplitude. What is left after the color 
bracketing is the spinor line connecting the two external fermions tt with all the possible ways and 
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the rest are the tensor integrals and scalar products of external momenta produced by propagators 
outside the loop and finally the polarization vectors of the gluons. 

The spinor line contribution is calculated by Fortran. After manipulating it using all the pos- 
sible contracting equations and simplifying the result by moving the appropriate momenta to the 
corners of the spinors and applying the Dirac equation, we end up with the simplest structure. Then 
we glue the two spinors left and right and calculate the result numerically in Fortran. 

The tensor integrals are identified and flagged by Form. There is a Fortran program that cal- 



culates them according to the New Recursive Reduction Method of [12]. The program called 
OLOTIC (One LOop Tensor Integral Calculator), applies the formulas of the paper expressing the 
tensor integrals in terms of scalar master integrals. In a few words, the idea of the reduction is 




IJ,stuv 



Figure 1: Graphical presentation of the tensor reduction procedure 

that, one-loop n-point rank R tensor Feynman integrals [in short: («,/?) -integrals] for « < 6 with 
/? < n are represented in terms of {n,R— 1)- and {n—l,R— l)-integrals (see e.g. Figure [l]). The re- 
cursive scheme is very convenient for explicit calculations and contains the complete calculational 
chain of tensor reduction, for both massive and massless propagators, and works with dimensional 
regularization. The general case of tensor integrals using dimensional regularization was treated 
in a series of papers [13, [I^, |^, |l6[ 17], thereby allowing also for massless particles. The only 



necessary package to be added is one for the evaluation of 1 -point to 4-point scalar integrals. These 
are calculated by QcdLoop [[HI]. This is a very interesting property of the reduction, because us- 
ing Fortran and switching off the relative scalar master integrals we have all the information that 
is needed to determine the coefficients of the scalar master integrals that arise in calculations by 
unitarity methods. In the same way, by switching off all the scalar master integrals we can also 
extract the rational term. Finally, due to the capability of QcdLoop to give the l/e and l/e^ terms, 
we can cross check the correctness of our loop amplitude as these terms are proportional to the tree 
level amphtude. 

The stripped amplitude in total is optimized by Maple. We divide the total number of diagrams 
to groups of hexagons, pentagons, boxes, triangles, and bubbles, and then we translate the output 
to Maple. Maple optimizes it and feeds it to a Fortran subroutine calculating the amplitude contri- 
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butions of the above groups. As mentioned above, these structures are single colored ones. So each 
one of them is multiplied with the appropriate color structure stored in the color subroutine. 

Connecting all the pieces together, the stripped amplitudes including tensor integrals and 
spinor lines separately calculated in specific subroutines, the color structures precalculated and 
stored in an array, we calculate the total loop amplitude. In the same way, using much simpler 
routines, we calculate the tree level case. 



DIANA 

Diagram 
construction 
Output (form) 



MAPLE INPUT 

cRanl<0.m(1...4) 

cRank5.m(1...4) 



gm{line,n1,...,n9) 
Spinor structures 



MADGRAPH 
momenta.dat 



fmmm 



bub m.frm 



Color2fortran.frm 
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FORTRAN OPT 

cFi_rtSum3(1...4) 

cS_rtSum23(1...4) 



IVIain fortran program 



50 different 

Structures 

color. F 




Hex(Sum_6(4)) 
Bubble(Sum_2(4)) 



Figure 2: Graphical presentation of the case of gg -^ ttgg 



The biggest advantage of the program is the fact that it is flexible. With minor changes, we can 
implement it for calculating different processes. We just have to provide Diana with the process 
that we are interested in, fix the Feynman rules and be careful with manipulation of possible spinor 
lines and colors. All the rest, tensor reduction and optimization of the code are done automatically 
by the program. 
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